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\ Abstract. Measurements of jet shapes and subjet multiplicities in photoproduction 

performed by ZEUS are presented and compared to leading-logarithm parton-shower 
Monte Carlo models. The predicted differences on the size of gluon- and quark-initiated 
jets are used to select samples to study the dynamics of the subprocesses. 
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§ : INTRODUCTION 

The internal structure of a jet depends mainly on the type of primary parton 
(quark or gluon) from which it originated and to a lesser extent on the particular 
<j_) ; hard scattering process. QCD predicts that (a) at sufficiently high transverse energy 

of the jet (E^), where fragmentation effects become negligible, the jet structure is 
driven by gluon emission from the primary parton; and (b) gluon jets are broader 
than quark jets due to the larger colour charge of the gluon. In the first part of 
this article, the measurements performed by ZEUS [1] to test the QCD prediction 
^ ■ (a) are presented and compared to leading-logarithm parton-shower Monte Carlo 

models. The lowest non-trivial order contribution to the jet substructure is given 
by 0(aa 2 ) calculations and, therefore measurements of jet substructure provide a 
stringest test of QCD beyond leading order (LO). The results of taking prediction 
(b) into account to test the dynamics of the subprocesses [2] are shown in the 
second part of this article. 

At HERA, positrons of energy E e = 27.5 GeV collide with protons of energy 
E p = 820 GeV. The main source of jets is hard scattering in 7^ interactions in 
which a quasi- real photon (Q 2 « 0, where Q 2 is the virtuality of the photon) 
emitted by the positron beam interacts with a parton from the proton to produce 
two jets in the final state. In LO QCD, there are two processes which contribute 
to the jet production cross section: the resolved process (figure la) in which the 
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(a) (b) 
FIGURE 1. (a) Resolved- and (b) direct-photon processes. 

photon interacts through its partonic content, and the direct process (figure lb) in 
which the photon interacts as a point-like particle. 

In the kinematic regime studied here, the dominant subprocesses are q 1 g p — > qg 
in resolved and 75 — > qq in direct. The kinematics of these processes are such that 
the majority of the jets in the region of pseudorapitity-f (rf^) below originate 
from quarks and the fraction of gluon-initiated jets increases as rf ei increases. 



JET SUBSTRUCTURE 



The internal structure of jets can be studied by means of the jet shape and the 
sub jet multiplicity. 



Jet shape 

The differential jet shape is defined as the average fraction of the jet's transverse 
energy that lies inside an annulus in the pseudorapidity (77) - azimuth (tp) plane 

of radii r ± Ar/2, where r = ^ (Arf) 2 + (A<^) 2 , concentric with the jet axis and is 
given by 

n J \_ v ffr(r-f ,r + f ) 

P[r) N jets 'Ar'^ s E T (0,r = l) ' 

where Et(t — Ar/2,r + Ar/2) is the transverse energy within the given annulus 
and Nj ets is the total number of jets in the sample. The integrated jet shape is the 
average fraction of the jet's transverse energy that lies inside a cone in the r] — <p 
plane of radius r concentric with the jet axis, 



% The pseudorapidity is defined as rj = — ln(tan |), where the polar angle 9 is taken with respect 
to the proton beam direction. The ZEUS coordinate system is defined as right-handed with the 
Z-axis pointing in the proton beam direction, hereafter referred to as forward, and the X-axis 
horizontal, pointing towards the centre of HERA. 
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(a) (b) 

FIGURE 2. (a) Measured differential jet shape p(r) vs. r (black dots), (b) Measured integrated 
jet shape ip(r = 0.3) vs. rf et (black dots). PYTHIA Monte Carlo calculations are shown for 
comparison. 
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Measurements of the differential jet shape p(r) have been performed by ZEUS [1] 
using an inclusive sample of jets. The jets have been identified using the longitudi- 
nally invariant kx cluster algorithm [3] in the inclusive mode. The jets were required 
to have E^ 1 > 17 GeV and —1 < rf et < 2. The measurements were performed 
in the kinematic region 0.2 < y < 0.85, where y is the inelasticity variable, and 
Q 2 < 1 GeV 2 . Figure 2a shows p(r) as a function of r in different regions of rf et . 
The data show that the jets become broader as if et increases. 

The solid histograms in figure 2a are calculations from the leading-logarithm 
parton-shower Monte Carlo PYTHIA [4]. The predictions, which include initial 
and final state QCD radiation, give a good description of the data. The comparison 
of the predictions for gluon- (dashed histograms) and quark-initiated (dot-dashed 
histograms) jets to the data shows that the measured jets are quark-like for — 1 < 
r]^ et < and become increasingly more gluon-like as nf et increases. 

The quark and gluon content of the final state has been investigated in more 
detail by studying the rf et dependence of the integrated jet shape at a fixed value 
of r [1]. Figure 2b shows ip(r) for r = 0.3 as a function of nf et . The measured 
jet shape decreases with rf ei , i.e. the jets become broader as rf et increases. The 



comparison between the predictions of PYTHIA for gluon- and quark-initiated jets 
and the data shows that the broadening of the jets is consistent with an increasing 
fraction of gluon-initiated jets as rf et increases. 

Subjet multiplicity 

The application of the k T cluster algorithm allows the study of the internal 
structure of the jets by means of jet-like objects (subjets) for which resummed 
QCD calculations are possible. Subjets are resolved within a jet by reapplying the 
kr cluster algorithm until for every pair of particles the quantity 

dij = mm(E Ti , E Tj ) 2 ■ [fa - r]j? + (<A - Vjf] 

is above y cut ■ (E^) 2 . The theoretical advantages of this observable are that safe 
observables can be defined at any order in perturbative QCD and resummed calcula- 
tions are possible. Subjets are also a useful tool to investigate the colour dynamics. 
Finally, the uncertainties coming from hadronisation corrections are small when 
y cu t is not too small. 

The mean subjet multiplicity < n su bj e t > has been measured [1] for an inclusive 
sample of jets with E^ 1 > 15 GeV and —1 < rf et < 2 in the kinematic region given 
by 0.2 < y < 0.85 and Q 2 < 1 GeV 2 . Figure 3a shows < n su bj e t > as a function 
of ycut- The predictions of PYTHIA and HERWIG [5], which include initial and 
final state QCD radiation, provide a good description of the data. The measured 
< n su bjet > for a fixed y cut value of 0.01 (see figure 3b) increases as rf et increases. 
The comparison to the predictions for gluon- and quark-initiated jets shows that 
the increase in < n su bj e t > as nf et increases is consistent with the predicted increase 
in the fraction of gluon-initiated jets. This result is consistent with the one obtained 
from the measurements of jet shapes. 

SUBSTRUCTURE DEPENDENCE OF DUET CROSS 

SECTIONS 

The predictions of the Monte Carlo models for the jet shape and subjet multi- 
plicity reproduce the data well and show the expected differences for quark- and 
gluon-initiated jets. Therefore, the Monte Carlo events have been used to devise a 
method to select samples enriched in quark- and gluon-initiated jets. The samples 
are selected by exploiting the QCD prediction that gluon-initiated jets should be 
broader than quark-initiated jets. 

The predicted distributions of i])(r) for r = 0.3 and < n su bj e t > for y cut = 0.0005 
for quark- and gluon-initiated jets are clearly different (see figure 4). A sample 
enriched in quark-initiated ("thin") jets has been selected by requiring ip(r = 0.3) > 
0.8 and a sample enriched in gluon-initiated ("thick") jets has been selected by 
requiring ip(r = 0.3) < 0.6. Alternative samples to support the results obtained 
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FIGURE 3. (a) Measured mean subjet multiplicity < n su bjet > vs. y cu t (black dots), (b) 
Measured mean subjet multiplicity < n su bj e t(y C ut = 0.01) > vs. rf et (black dots). PYTHIA and 
HERWIG Monte Carlo calculations are shown for comparison. 

with the jet shape selection have been selected using the subjet multiplicity. Jets 
are called "thin" if n su bj e t(y C ut = 0.0005) < 4 and the jets are called "thick" if 
n su bjet(ycut = 0.0005) > 6. 

Measurements of da / drf et 

Using the jet shape selection into "thin" and "thick" jet samples, dijet cross 
sections have been measured [2] as a function of r/- 76 ' for events with at least two 
jets of E^ 1 > 14 GeV and —1 < rf et < 2.5 in the kinematic region given by 
0.2 < y < 0.85 and Q 2 < 1 GeV 2 . Figure 5 shows da/d^ et for "thick" and "thin" 
jets. The cross section for "thick" jets displays a very different shape than that 
of the "thin" jet sample: the rf et distribution for the "thin" jet sample peaks at 
about 0.7, whereas the "thick" jet sample distribution peaks at r]^ et ps 1.5. 

The predictions of PYTHIA and HERWIG (see figure 5a) give a good description 
of the data. The Monte Carlo distributions have been normalised to the total 
measured cross section of each type after applying the same jet shape selection as 
in the data. PYTHIA and HERWIG predict a similar parton content of the final 
state: the "thick" jet sample is composed of 15 — 17% of gg subprocesses in the final 
state, 54 — 58% of gq and 25 — 31% of qq. The "thin" jet sample contains 54 — 56% 
of qq, 41% of qg and 3 — 5% of gg. Therefore, the "thin" jet sample is indeed 
dominated by quark-initiated jets in the final state and the "thick" jet sample has 
a high content of gluon-initiated jets coming mainly from the final-state gluon of 
the subprocess q 7 g p — > qg. The measurements are compared to the predictions of 




FIGURE 4. The predicted integrated jet shape V( r = 0.3) (a) and subjet multiplic- 
ity n su bjet(ycut — 0.0005) (b) distributions at the hadron level for samples of gluon- and 
quark-initiated jets simulated using the program PYTHIA. 

quark- and gluon-initiated jets with the same jet shape selection as for the data 
in figure 5b. A "thin" -quark jet sample gives a good description of the "thin" 
jet sample, whereas the "thin" -gluon jet sample peaks one unit of pseudorapidity 
higher. A sample of "thick" -quark jets alone cannot describe the "thick"-jet sample 
and a large content of gluon-initiated jets is needed to reproduce the shape of the 
data. 

Measurements of da/d cos 6* 

The underlying parton dynamics is reflected in the distribution of the scattering 
angle in the dijet centre-of-mass system, 9*. The cos9* distribution is sensitive 
to the spin of the exchanged particle: for gluon exchange, the cross section is 
proportional to (1 — |cos#*|)~ 2 , whereas for quark exchange the cross section is 
proportional to (1 — | cos9*\)~ 1 . In a previous publication [6], the | cos#*| distribu- 
tion for resolved processes (which are dominated by subprocesses mediated by gluon 
exchange) was observed to rise more steeply than that of direct processes (which 
proceed via quark exchange) at high | cos 9* | values. In this analysis [2] , samples 
of jets are selected according to their internal structure. This method provides a 
handle to better understand the dynamics of the subprocesses. 

For samples of two "thick" or two "thin" jets, only the absolute value of cos#* 
can be determined because the outgoing jets are indistinguishable. The dijet cross 
section as a function of | cos#*| for samples of two "thick"-jet events and two "thin"- 
jet events has been measured for dijet invariant masses M JJ > 47 GeV. The cross 
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FIGURE 5. Measured da/drf> et for samples of "thick" (black dots) and "thin" (open circles) 
jets selected according to their shape. The thick error bars represent the statistical errors of the 
data, and the thin error bars show the statistical errors and uncorrelated systematic uncertainties 
added in quadrature. The shaded band displays the uncertainty due to the absolute energy scale 
of the jets. In (a) Monte Carlo calculations using PYTHIA and HERWIG for resolved- plus 
direct-photon processes and in (b) the predictions of PYTHIA for quark- and gluon-initiated jets 
are shown for comparison. 
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FIGURE 6. Measured da / 'd\cos9*\ for two "thick"-jet events (black dots) and two "thin"-jet 
events (open circles) selected according to their shape (a) or to their subjet multiplicity (b). 
PYTHIA and HERWIG Monte Carlo calculations are shown for comparison. 



sections are presented in figure 6 and are normalised so as to have a value of unity at 
low | cos 6**| to study the difference in slope as | cos#*| — > I. The | cos9*\ distribution 
for the two samples of dijet events increases as |cos#*| increases, however they 
exhibit a different slope. For comparison, PYTHIA and HERWIG Monte Carlo 
predictions are shown in figure 6. The predictions have been also normalised so 
as to have a value of unity at low | cos 9* | and give a reasonable description of the 
shape of the data. The different slope observed in the data can be understood in 
terms of the dominant subprocesses in the two samples: the two "thick" -jet sample 
is dominated by subprocesses mediated by gluon exchange (gg — > gg and qg — > qg) 
whereas the two "thin"-jet sample is dominated by subprocesses mediated by quark 
exchange {jg — > qq). 

The sample of events with one "thick" jet and one "thin" jet allows a measure- 
ment of the unfolded da/dcos6* cross section, since in this case the jets can be 
distinguished. Figure 7 shows the measured dijet cross section as a function of 
cos9*. The dijet angular distribution was measured with respect to the "thick" 
jet and shows a different behaviour on the negative and positive sides: the mea- 
sured cross section at cos9* = 0.7 is more than two times larger than the one at 
cos 9* = —0.7. The Monte Carlo models give a reasonable description of the shape 
of the measured da jd cos 9* (see figure 7). The measurements and the Monte Carlo 
calculations in figure 7 are normalised so as to have a value of unity at central values 
of cos 9* to study the difference in slope as cos 9* — > ±1. The observed asymmetry 
is understood in terms of the dominant subprocess: q^g p — > qg. The positive side is 
dominated by t— channel gluon exchange, whereas the negative side is dominated 
by u— channel quark exchange. 

SUMMARY AND CONCLUSIONS 

Measurements of jet shape and subjet multiplicity have been performed for in- 
clusive jet samples with E 3 ^ 1 > 15 GeV and — 1 < rf et < 2 in the kinematic regime 
given by 0.2 < y < 0.85 and Q 2 < 1 GeV 2 . The measured jet shape broadens and 
the mean subjet multiplicity increases as rj 3et increases. Leading-logarithm parton- 
shower Monte Carlo models with initial and final state QCD radiation give a good 
description of the data. The observed broadening of the jet shape and the increase 
in the mean subjet multiplicity as i] jet increases is consistent with an increase of 
the fraction of gluon jets. 

The Monte Carlo models reproduce the measurements of the jet shape and sub- 
jet multiplicity and display the expected differences for quark- and gluon-initiated 
jets and allow the use of the jet shape and subjet multiplicity to select samples 
enriched in quark- and gluon-initiated jets to study the dynamics of the subpro- 
cesses. Measurements of the dijet cross section as a function of rf et for samples of 
"thick" and "thin" jets show the expected behaviour for samples enriched in gluon- 
and quark-initiated jets. The | cos#*| distribution for a two "thick"-jet sample dis- 
plays a similar behaviour to the one expected for a sample enriched in processes 
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FIGURE 7. Measured da/dcos9* for one "thick" jet and one "thin" jet in an event (black dots) 
selected according to their shape (a) or to their subjet multiplicity (b). PYTHIA and HERWIG 
Monte Carlo calculations are shown for comparison. 



mediated by gluon exchange, whereas that for a two "thin"-jet sample shows a 
behaviour similar to a sample enriched in processes mediated by quark exchange. 
Finally, the cos#* distribution for a sample of events with one "thick" and one 
"thin" jet exhibits a large asymmetry consistent with the expected dominance of 
t— channel gluon (u— channel quark) exchange as cos#* — > +1(— 1). 
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